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Abstract. - We investigate the shape of the spectrum and the spectral fluctuations of the 
fe-body Embedded Gaussian Ensemble for Bosons in the dense limit, where the number of 
Bosons TO, — > oo while both k, the rank of the interaction, and I, the number of single-particle 
states, are kept fixed. We show that the relative fluctuations of the low spectral moments do 
not vanish in this limit, proving that the ensemble is non-ergodic. Numerical simulations yield 

h^ , spectra which display a strong tendency towards picket-fence type. The wave functions also 

f^ ■ deviate from canonical random-matrix behaviour. 



Introduction. - Random-matrix theory (RMT) successfully describes the statistical be- 
haviour of spectra and wave functions of a large variety of systems such as atoms, molecules, 
atomic nuclei and quantum dots [0,0]. However, this RMT modeling is not completely re- 
alistic since all many-body systems are effectively governed by one- and two-body forces. 

Q , This fact led to work on the two-body random ensembles for Fermions ||§-^ and to the in- 
troduction of the A:-body embedded ensembles by Mon and French 0. In the embedded 
ensembles, many-body states are constructed by distributing m particles over I degenerate 
single-particle levels. The matrix of the fc-body interaction with fc < m is taken in this basis. 

3 , For k < m, the 77i-body matrix elements of the random fc-body interaction are correlated: 
The number of independent random variables is smaller than in RMT. Do these more realistic 
embedded ensembles yield the same results as RMT? Early numerical simulations for inter- 
acting Fermion systems P,0] of rather small matrix dimension have shown that the spectral 
fluctuation properties of the embedded ensembles agree with those of RMT. Similar results 
were obtained in numerical simulations for Bosonic systems |8|^. Moreover, the Fermionic 
ensembles were shown to be ergodic. As far as we know, there are no results on the spectral 
ergodicity for Bosons. 

Recently, three of the present authors introduced a novel analytical approach to the 
Fermionic embedded ensembles in the limit of infinite matrix dimension (Z — > oo) |Ig|. The 
main results of this approach are: (i) For 2fc > m, the average spectrum has semicircle shape, 
and the spectral fluctuation properties coincide with those of RMT; (ii) the spectral density 
changes shape at or near 2k = m and becomes Gaussian; (iii) in the dilute limit {k <^m <^l) 
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the spectral fluctuations are completely uncorrelated (Poissonian); (iv) the spectral fluctua- 
tions change gradually from Wigner-Dyson to Poisson. 

In this Letter, we extended our work to the case of Bosons. As in the Fermionic case, 
we consider the limit of infinite matrix dimension. For Bosons this limit is realized either by 
letting I — > oo (the same limit as for Fermions), or by letting tti — > cx) with k and I fixed. 
This second case, the dense limit, is novel and has no analogue in the Fermionic case. In 
the following, we focus attention exclusively on the dense limit. We prove analytically that 
in this limit, the ensemble is not ergodic. Numerical results for the spectral correlations are 
obtained by both ensemble unfolding and by spectral unfolding. In the latter case we find 
a strong tendency of the spectra towards picket-fence behaviour. We also show that some 
eigenfunctions display Fock-space localization. Details of the derivations and a complete 
treatment including the limit Z — > oo, are given in Ref. |ll| . 

Definitions. - We consider m spinless Bosons in I degenerate single-particle states with 
associated creation and annihilation operators b- and bj where j = l,...,l. Hilbert space 

is spanned by the N — (^ ^™^ ) orthonormal m-particle states |/i) {fi — 1,. . . ,N), written 

as l/x) = [M{ji, . . . , jm)]"^^]i ■ ■ ■ ^L ^'^^^ ji < 32 < ■■■ < im and M{ji, .■.,3m) a normal- 
ization constant. The states |/i) are equivalently characterized by a sequence of occupation 
numbers (ni, . . . , n;) of the I single-particle states with &I ■ • • b*- = (fe{)"^ . . . (&/)"' ■ Setting 
Af{ji, . ■ . ,jm) = V"i! ■ . .nil normalizes all ?Ti-particle states \fj) to 1. 

The bosonic m~particle states are coupled through a random fc-body interaction Vk{P) 
with k < m, given by 

■5-^ 61 . . .b\ bi. . . . bu 

'-''^^ ^ ■■'■ '■' -mtrminr^)^ « 

We refer to k as to the rank of the interaction. As in the canonical case, the labels /3 — 1 
and (5 = 2 denote the orthogonal and the unitary ensemble, respectively. The matrix element 
'^h,---,3k\ii,---,ik of the fc-body interaction taken between the single-particle states j'l, . ■ . ,jk and 
«i, . . . ,i/c is totally symmetric with respect to ji, . . . ,jk and ii, . . . ,ik. The elements differing 
in the sequence of indices {j'l . . . jk]ii . . . ik} (except for permutations of {j'l . . . jk} and of 
{ii . . .ik} and for symmetries specified by /3) are uncorrelated Gaussian-distributed random 
variables with zero mean and a common second moment Vq . Without loss of generality we put 
iiq = 1 in the sequel. The normalization coefficients A/" (j'l, ... ,jfc) 'And Af{ii, ... ,ik) in Eq. (|l|) 
guarantee that for k = m, the embedded ensembles are identical to the canonical ensembles 
of random-matrix theory. This defines the Bosonic fc-body embedded Gaussian orthogonal 
(unitary) ensemble of random matrices, respectively, in short BEGOE(fc) and BEGUE(fc). 
The number of independent random variables is given by K/j = /?( ^j,~^) [( ^ k~^) ~^ <^/3i]/2- 

The Second Moment. - By virtue of the randomness of 14- (/3), the elements of the 
matrix {i>\Vk{P)\fJ.) are random variables with a Gaussian probability distribution and zero 
mean value. The spectral properties are completely determined by the second moment 

Bi%AP) = TiAMWiWkUm 

a(k),'y{k) 

The overbar denotes the average over the ensemble. We have simplified the notation by 
introducing the operators B*,. = [Af{ji, . . . ,jr)]~^b', . . .b', and the adjoints Ba(r)- The 
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index a{r) is a short-hand notation for the rank r and for the sequence of indices {j'l . . . jV}- 
The matrix B^^'{P) is Hermitean in the pairs of indices (/i, i^) and {p,(j). It is easy to prove 
the "duahty" relation Bl,J^pa{'2) = Bj^^^' (2) which for /? = 2 connects the second moments 
of the fc-body and the {m — fc)-body interactions. For k = m and (3 = 2, the right-hand 
side of Eq. (||) reduces to S^u^pa, and correspondingly for /3 = 1. This shows that for k = m, 
BEGUE(fc) and BEGOE(fc) reduce to the GUE and GOE, respectively. 

In the Fermionic case, results for the shape of the spectrum and for the spectral fluctuations 
were obtained with the help of the eigenvector decomposition of the second moment. A similar 
decomposition exists in the case of Bosons and is derived in Ref. pi. In the following, we only 
state the results of this derivation. We solve the eigenvalue equation '}2na^\^'^-p<^i'^)^^p — 
A^^^{k)CpT , with eigenvectors C^*°^ and with eigenvalues A(*^(fc) given by 

A(s)(fc)^(^--^^(^^ + - + -^-l^ (3) 

Here, s = 0, . . . ,m while a labels the degenerate eigenvectors. For the degree of degeneracy 
£)(*') of the eigenvalues A(")(A:) we find D^o) = 1 and D^"') = ('+7^)^ - Ct-7^)^ for s > 1- It 
follows that X^rio -^''''' ~ ^^ showing that the eigenvectors form a complete set. We choose 
Hermitean linear combinations of the degenerate eigenvectors which obey the orthonormality 
condition Xluj/ C'/jtf W/i ' = NSstSab- These results can be extended to f3 = I. Hence, the 
matrix B^^\f3) can be expanded as 

BjS,pAP) = ^ E ^^'^(^) Et^^t^^^^r^ + spiCifciT'^] . (4) 

Low Moments ofVk- - Using the eigenvector decomposition of B^^\(S), duality and the 
orthonormality of the C^*"-''s, explicit expressions for the low moments of Vfe(/3) are obtained. 
From these we calculate three ratios that yield information on the shape of the spectral density. 
The ratio 5* measures the fluctuations of the center of the spectrum in units of the average 
width of the spectrum. The ratio i? measures the relative fluctuations of the width of the 
spectrum. The ratio Q is related to the kurtosis k = Q + 2 and marks the difference between 
the semicircular (Q = 0) and the Gaussian shape (Q = 1). The definitions of S*, i?, and Q are 
given in Ref. | [lC| ]. 

We are particularly interested in the behaviour of the three ratios in the dense limit, i.e. 
the case iV — > cx) attained by letting rn — > oo while keeping both / and k fixed. We obtain 

hm 5(fc,m,o = — (i+^^i)a'')rr^)'' ^^^ 

lim R{k,m,l) = ^s=ol\.k+JK k+s I J L W n ^ ^g^ 

[(fc)+^/5iE.=o(fe+J( k+s ) d('' 

with S'^^ = ( ^s^^) • I^'l- (0) applies only in the unitary case (/3 = 2) and implies that in the 
dense limit, the average spectrum has Gaussian shape. We have not been able to extend this 
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Fig. 1 - Spectral density of the BEG0E(2) for I — 2 and m = 3000, normalized to the dimension 
A'^ = 3001 of Hilbert space. The curves (a) and (b) show the level densities of two members of 
the ensemble; curve (c) shows the ensemble-averaged spectral density. The dashed curve shows the 
theoretically predicted Gaussian shape of the average spectrum. 

Fig. 2 - Nearest-neighbour spacing distribution P{s) obtained by ensemble unfolding (filled his- 
togram) and spectral unfolding (empty histogram). The dotted curve gives the Wigner surmise and 
the dashed one the Poisson distribution. 



analytical result to the orthogonal case (3=1. On physical grounds, however, an equation 
analogous to Eq. (0) is expected to be valid, also for the BEGOE(fc). We conclude that the 
shape of the average spectrum has Gaussian form. This is in keeping with the results of 
Ref. |l|. 

A more important and surprising result lies in the fact that the right-hand sides of Eqs. (^ 
and (g) do not vanish: The fluctuations of the centroids and of the variances of individual 
spectra do not vanish asymptotically. This feature differs from the behaviour both of canonical 
RMT and of the embedded Fermionic ensembles. We are led to the important conclusion that 
the Bosonic ensembles are not ergodic in the dense limit m ^ oo with k and I fixed. Non- 
ergodicity appears to be a consequence of the fact that the number of independent random 
variables Kp in the ensemble does not grow with m, but stays finite for fixed k and I. 

Numerical results. - For lack of analytical techniques, we use numerical simulations to 
obtain infromation on the spectral fluctuation properties of BEGOE(fc) in the dense limit. 
With I = 2 and k = 2, the dimension of Hilbert space is given by A^ = m -I- 1, and m, ^ Hs 
easily attainable numerically. The number of independent random variables is Ki = 6. The 
many-particle states |/i„) are written as (m — n,n) where m — n and n indicate the number 
of Bosons occupying the first and the second single-particle state, respectively. In the basis 
{|/io), |mi)i • ■ • I l/^m)}, the Hamiltonian matrix attains band structure, with non-zero matrix 
elements on the main diagonal and on the k closest diagonals. 

Fig. |l| shows the level densities of two members of the ensemble and the ensemble-averaged 
density obtained from 1512 spectra {m = 3000). We have also plotted a Gaussian density 
whose width is given by the theoretical prediction [Q and which agrees with the shape of the 
ensemble-averaged density. The striking differences between each of the two spectral densities 
and between these and the ensemble-averaged spectral density illustrate the non-ergodic 
character of BEG0E(2) for m ^ Z. The Gaussian form of the ensemble-averaged spectrum 
arises as a consequence of the Central Limit Theorem. It has no bearing on individual spectral 
shapes. 
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Fig. 3 - As-statistics measured at the centers of the spectra after ensemble unfolding (thin solid line) 
and spectral unfolding (bold line). For comparison we have plotted the results for the GOE (dotted 
curve), for a Poissonian spectrum (dashed line), and for a picket-fence spectrum (dotted-dashed line). 

Fig. 4 - Detail of the staircase function rj{E) for a typical member of BEGOE(2) with m = 1000 and 
I = 2. The vertical lines help to show how two (almost) equidistant spectra overlap around the level 
620. We note the change in the density of states after this level. 



The evaluation of measures of spectral fluctuations requires the spectra to be unfolded. 
Unfolding can be done either individually for each spectrum (spectral unfolding) or by a 
single transformation which is applied to all spectra of the ensemble (ensemble unfolding). 
We employed both types of unfolding in our statistical analysis. For spectral unfolding, a 
polynomial was fitted to the staircase function of each realization of a spectrum. Starting 
from the value unity, we changed the degree of the polynomial for each realization until the 
first minimum of the associated y^ distribution was reached. The maximum degree considered 
was 20. The first minimum was typically found around degree 11. Ensemble unfolding was 
carried out by first averaging the staircase functions over the ensemble. The averaged staircase 
function was fitted by a polynomial of degree 11. We have analyzed 1512 realizations for 
m = 3000. After fitting the polynomials by including all levels of each spectrum for both 
spectral and ensemble unfolding, we only considered 60% of all levels, i.e., 1800 levels for each 
realization located in the centre part of each spectrum. 

In Figs, g and || we show the nearest-neighbour spacing distribution P{s) and the A3- 
statistics obtained for m = 3000 by ensemble unfolding and by spectral unfolding. In both 
cases the distribution P{s) corresponds neither to the Wigner surmise nor to a Poisson dis- 
tribution. In the case of ensemble unfolding the level repulsion characteristic of the GOE is 
clearly lost. For spectral unfolding the distribution P{s) is dominated by a prominent peak 
centered at s = 1 suggesting that individual spectra have an almost constant level spacing. 
The As-statistics obtained by ensemble unfolding deviates rapidly from GOE behaviour and 
increases almost linearly. In the case of spectral unfolding, A3(L) is almost constant up to 
L ~ 20. This result again suggests a tendency of individual spectra towards a picket-fence-like 
behaviour. Beyond this range, /^^{L) grows albeit less rapidly than for ensemble unfolding. 

We turn to the detailed structure of individual spectra. Fig. ^ shows the staircase function 
'ri{E) for a typical member of the ensemble {k = 2, I = 2, m = 1000). The spectrum is 
dominated by levels with almost constant spacing. The staircase function typically displays 
one or more points where an abrupt change in the density of states takes place. In Fig. 0, the 
spectrum has (almost) constant spacing up to the level 619. From level 620 on, the spacing 
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Fig. 5 - Probabilities \cinf for eigenvectors belonging to eigenvalues in the vicinity of the kink shown 
in Fig. W. The two overlapping segments of nearly equidistant levels are easily distinguished by the 
structure of the eigenfunctions. 



between neighbouring levels is no longer constant. However, the spacing between next-to- 
nearest-neighbours does have this property and is almost the same as the nearest-neighbour 
spacing before level 620. This is illustrated by the vertical lines plotted above and below 
the staircase function. These observations imply non-stationary properties of the spectra. 
In addition, our result suggests that the spectrum of an individual member of the ensemble 
consists of pieces of overlapping segments of spectra with almost constant level spacings. The 
kink at level 620 in Fig. |j marks the left edge of the overlap region. 

The overlap of two segments of almost equidistant levels also influences the structure of the 
eigenfunctions. We write the eigenvectors \i) as linear combinations of the ordered many-body 
basis states |/i„), |i) = J2n=o CmiMn)- In Fig- H we plot the probabilities |ci„p of eigenvectors 
\i) belonging to eigenvalues in the vicinity of the kink shown in Fig. H. Eigenvectors up 
to i = 619 behave similarly and are somewhat extended, although the distribution of the 
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intensities clearly deviates from the Porter-Thomas distribution expected from RMT. This 
behaviour changes abruptly aii — 620. We emphasize the difference in the scales used in Fig. 
for i = 619 and for the even i values starting with i = 620. The eigenvectors with i = 621, 623 
display the same behaviour as the eigenvectors up to i = 619 and, thus, correspond to the 
first segment of the two equidistant spectra. The eigenvectors with i — 620, 622, 624 are 
much more localized in Fock space. They differ in the number of intensity oscillations. As 
i increases, so does this number, and the levels on this second segment become more and 
more delocalized. At some point the spread of the eigenfunctions in the second segment is 
indistinguishable from that in the first one. 

Conclusions. - We have studied the spectral behaviour of the Bosonic embedded ensem- 
bles in the dense limit attained by letting m — > cxd with I and k fixed. We have shown that in 
this limit, both the BEGOE(A;) and the BEGUE(fc) are non-ergodic. To the best of our knowl- 
edge, this is the first case of a random-matrix model for which such non-ergodic behaviour 
has been established in the limit of infinite matrix dimension. Moreover, we have shown that 
the spectral fluctuations deviate strongly from RMT results. This result disagrees with the 
conclusions based on numerical simulations of Refs. |§|,||]. We ascribe this disagreement to 
the fact that the ratios m/l considered in Refs. [g,^ were too small to produce deviations 
from Wigner-Dyson statistics. Our numerical simulations provide evidence for both a non- 
stationary and a picket-fence-type of behaviour in both short- and long-range correlations 
of the spectra. More precisely, individual spectra seem generically to consist of overlapping 
segments of spectra of picket-fence type. In the overlap region, the eigenfunctions of the lev- 
els in the two segments differ markedly, the eigenfunctions in the second segment displaying 
strong localization in Fock space. We have not yet attained an analytical understanding of 
these properties which we attribute to the small number of independent random variables Kfj 
(with Kfs independent of m) , and to the specific structure of the Hamiltonian matrix of the 
Bosonic ensembles. In particular, we do not know how the spectra change as I increases. 

* * * 
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